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We study quantum states of electrons in magnetically doped quantum dots as a function of ex- 
change coupling between electron and impurity spins, the strength of Coulomb interaction, confining 
potential, and the number of electrons. The magnetic phase diagram of quantum dots, doped with 
a large number of magnetic Mn impurities, can be described by the energy gap in the spectrum 
of electrons and the mean field electron-Mn exchange coupling. A competition between these two 
parameters leads to a transition between spin-unpolarized and spin-polarized states, in the absence 
of applied magnetic field. Tuning the energy gap by electrostatic control of nonparabolicity of the 
confining potential can enable control of magnetization even at the fixed number of electrons. We 
illustrate our findings by directly comparing Mn-doped quantum dots with parabolic and Gaussian 
confining potential. 
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Recent experimental 
oretical 



II i, 0, i, i and the- 
^ ^ studies reveal that semi- 
conductor quantum dots (QDs) [3, [l3| are desirable for 
tailoring magnetism which persists at much high tem- 
peratures than in their bulk counterparts |18| . These 
nanoscale realizations of dilute magnetic semiconductors 
(DMS) may be suitable for a versatile control of 
spin and magnetism with potential applications in infor- 
mation storage technology and spintronics ^0]. An im- 
portant property of the carrier-mediated ferromagnetism 
in bulk-like DMS is the optical and electrical control of 
the ma gnet ic ordering by changing the number of carri- 



ers [21|, l22|, [23] ■ However, because of the strong inter- 



play between quantum confinement and the many-body 
electron-electron (e-e) Coulomb interactions in QDs, it 
may be possible to achieve and control the magnetic or- 
dering of carrier spin and magnetic impurities, with a 
fixed number of carriers 

In this paper, we present a theoretical analysis to de- 
scribe the quantum states of the electrons confined in 
II- VI quantum dots doped with Mn. In (II,Mn)VI ma- 
terials, Mn does not change the number of electrons (N) 
because it is isoelectronic with gronp-II elements. In 
QDs, however, additional carriers are controlled by either 
chemical doping or by an external electrostatic potential 
applied to the metallic gates. Even for a small number of 
interacting electrons, the study of QDs doped with large 
numbers of Mn (> 10) becomes computationally inacces- 
sible to exact diagaonalization techniques For 
this reason we employ here a mean field theory of mag- 
netic QDs in which Mn acts like an external magnetic 
field and the weak coupling between electrons and mag- 
netic impurities can be described as spin-spin exchange 
interaction [l^. Because electrons are confined in all 
three dimensions in nanometer-sized QDs, we anticipate 
strong electron-electron (e-e) Coulomb interaction. We 
use a finite temperature local spin density approxima- 
tion (LSDA) [iSl to incorporate e-e Coulomb interaction 
in this work. 



The wave function of electrons in QDs can be expanded 
in terms of its planar ipiaip) with p = {x,y), and sub- 
band wave function S,{z). The effective two-dimensional 
(2D) Kohn-Sham (KS) equations H^i„{p) = eiaipiaip), 
in LSDA can be obtained by integrating ^(z) 25|, assum- 
ing that the first subband is filled. Here e^cr are the KS 
eigenenergies, a — +1 and —1 for spin up (|) and down 
(I), and H is the Kohn-Sham Hamiltonian which can be 
expressed as 



Vqd + iVh + iV^c - ^h,d{p). (1) 



where h is the Planck constant, m* is the electron ef- 
fective mass, and Vqd = Vq exp{— / 5^) is the 2D 
Gaussian confining potential of the quantum dot ( ID 
parabolic confinement is chosen along the z-axis). In 
Eq.(IT]), Vh and V^c are Hartree and spin dependent 
exchange-correlation potentials [23| , while 7 accounts for 
reduction of Coulomb strength due to screening effects of 
the gate electrodes i27j and 



hsdip) = Jem / dz\^iz)\'^B 



'M 



( Mh{p,z) 



(2) 



Jem = Jscin-mM is the e-Mn exchange coupling, Jsd is the 
exchange coupling between electron and single Mn, Um is 
the spatial-averaged density of Mn, and M = 5/2 is the 
spin of Mn. Bm{x) is the Brillouin function [2^ in which 
the argument x contains kB the Boltzmann constant, T 
the absolute temperature, and 



J. 



b{p, z) = -J^i {Mdfi, z)) + ^K(p, z) - n^{p, z)], (3) 



the effective field seen by the Mn [15|. The first term 
in b{p, z) describes the mean field of the direct Mn-Mn 
antiferromagnetic coupling [l^ . with 



{M,{p,z)) = MBm 



( Mb{p,z) 



(4) 
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FIG. 1: The schematic single particle levels of a 2D parabolic 
(a) and 2D Gaussian (b) confining potentials corresponding 
to dashed and bold lines in (c). 



and 



Ap,z) = J2\i,M\'\az)\'f{e^.l 



(5) 



is the spin-resolved electron density, /(e) = l/{cxp[(e — 
fi)/kBT] + 1} is the Fermi function, and fi is the 
temperature-dependent chemical potential, and is the 
self-consistent solution of iV = ^ f[€-ia) which ensures 
that N — + Ni'is an integer number. 

To obtain solutions of Eqs. (P)-©, we consider 
(Cd,Mn)Te QD with Vb = -128 meV, and 5 = 15.9 nm, a 
perpendicular width of 1 nm, J^d = 0.015 eV nm'^, n,„ = 
0, 0.025, 0.1 nm~^, corresponding to the approximate 
number of Mn, Ny^ = 0, 45, 180, Je,„ = 0, 0.94, 3.75 
meV, and J^^ = 0, 0.005, 0.02 meV, respectively. For 
CdTe, considered here, a*g — 5.29 nm, and Ry* — 12.8 
meV are the effective Bohr radius and Rydberg energy, 
while the additional material parameters are to* = 0.106, 
and e = 10.6 

The shell structure of the 2D parabolic (a) and 2D 
Gaussian (b) potentials are shown in Fig. [TJ The cor- 
responding confining potentials are shown in Fig. [ijc) 
where the dashed and bold lines represent the 2D 
parabolic and 2D Gaussian potentials. The distinction 
between parabolic and Gaussian potentials can be ex- 
perimentally resolved by photo-emission spectroscopy or 
transport measurements of single particle states in ex- 
ternal magnetic field [ii[l3. The energy gap between 
S-, p-, and d-orbitals is characterized by loq. For a 2D 
Gaussian potential, the characteristic energy associated 
with the perpendicular confinement, cjq, is calculated by 
expanding Vqd in the vicinity of the minimum which 
yields Vqd — Vb -I- TO*Wgp^/2 + . . ., with the strength 
uoq = yj2\Vo\/m* /5. Because of the circular symmetry of 
the Gaussian potential, the z-component of the angular 
momentum Iz is a good quantum number. As a result 
of the geometrical symmetry of the confining potential, 
the states in p-shell with Iz — ±1 are degenerate. In con- 
trast to the 2D parabolic potential 3, 17 1, where the 
dynamical symmetry of the confining potential implies 
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FIG. 2: Single particle levels of the d-shell as a function of 
effective Zeeman energy hsd- The insets (a) and (b) schemat- 
ically show the ground state of A'^ = 8 depending on the rel- 
ative magnitude of hsd and the single particle energy gap A. 
At A = hsd (shown by dashed line) a transition from anti- 
ferromagnetic to ferromagnetic state is predicted. Filled and 
dashed lines represent the single particle eigenenergies of the 
electrons with spin up and down. 



occurrence of accidental degeneracies of the single parti- 
cle levels with different |^z|, the single particle levels in 
(i-shell are not completely degenerate. Degenerate levels 
and d_ are separated by an energy gap (1.5 meV) 
from the do-level, where the indices ±, refer to angular 
momentum Iz = ±1, 0. We note that the ordering of the 
single particle s-, p-, and d-eigenenergies in the Gaussian 
potential are similar to the ordering of the first few single 
particle eigenenergies in a square potential. 

In the absence of direct antiferromagnetic coupling be- 
tween Mn, and e-e Coulomb interaction, an identification 
of the magnetic states of the QDs based on the present 
mean field theory can be expressed in terms of two en- 
ergy scales: the single particle energy gap A, and the 
effective e-Mn exchange coupling Jem- The former is the 
energy difference between the highest occupied molecu- 
lar orbital (HOMO) and the lowest unoccupied molecular 
orbital (LUMO). It is a function of the number and spin 
of electrons, and the geometrical and dynamical symme- 
tries of the QD confining potential. Therefore A, which 
depends on the quantum confinement, can be tuned and 
controlled by the electric voltage applied to the metallic 
gates. On the other hand, Jerm the strength of the ef- 
fective Zeeman energy, hsd = {hsd{p)), which polarizes 
spin of electrons, is independent of the quantum con- 
finement. Thus the magnetic phase diagram of the QDs 
can be effectively described in terms of A, and Jem- In 
the absence of magnetic field we can then distinguish be- 
tween spin-polarized and spin-unpolarized states in QDs, 
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FIG. 3: Spin up (bold lines), and spin down (dashed lines) 
show Kohn-Sham eigenstates of N = 8, 9, 10 electrons in d- 
shell in 2D Gaussian confining potential with Vb = —128 meV, 
S — 15.9 nm, rim. = (a), and Um = 0.1 nm"'' (b). For clarity, 
s- and p- levels are not shown. 



conventionally also referred to as the ferromagnetic (FM) 
and antiferromagnetic ( AFM) states 13, l3i M \ ■ For ex- 
ample, the ground state of the electrons exhibits FM or- 
dering if Jem > A. The dependence of the single particle 
eigenenergies of the d-shell with iV = 8 on hgd is shown in 
Fig. [21 At hsd = the valence electrons with spin up and 
down occupy single particle level do, following the Pauli 
exclusion principle, and the ground state of = 8 ex- 



hibits AFM ordering with — — 4. With increasing 
the gap between spin up and spin down opens up. At 
hsd = A = 1.5 meV, where A = -Ed+Cd.) — Ed^, the first 
spin flip occurs. Because of well separated p- and d- or- 
bitals the second spin flip occurs at much higher energies 

as hsd = CJQ ~ 7 Jem. 

In LSDA, electrons are replaced by independent quasi- 
particles (electrons dressed by e-e interaction) , where KS 
eigenenergies are the energy levels of the quasi-particles. 
In this picture, e-e interaction rearrange the structure 
of quantized levels, and affects the shell structure of the 
QDs. In this limit, the quasi-particle gap A* (renor- 
malized A) is the energy difference between HOMO and 
LUMO of the KS eigenenergies. It is an appropriate pa- 
rameter, which can supersede the single particle gap A 
in identification of the magnetic states of the QDs. 

The KS levels in the d-shell of the 2D Gaussian poten- 
tial with fim = and A^ = 8, 9, 10 are shown in Fig. ^a). 
Because of circular symmetry of the confining potential, 
the energy levels of d+ and d- are degenerate. The en- 
ergy of do which was lower than d+ and d- due to the 
nonparabolicity of the Gaussian potential, increases by e- 
e Coulomb interaction, such that it becomes the highest 
energy level in d-shell, and A* changes sign. In contrast 
to non-interacting electrons, the ordering of the KS lev- 
els in S-, p-, and d-shells is closer to the one in circular 
potential than in the square potential. In the absence of 
magnetic impurities, as it is shown in Fig. [3Ja), N = 8 




FIG. 4: Energy of quantized states in 2D Gaussian confining 
potential with spin up (bold lines) and spin down (dashed 
lines) in d-shell as a function of effective Zeeman energy hsd 
corresponding to A* < (a), A* = (b), and A* > (c). 
Energy levels marked with circles are degenerate and shifted 
with respect to each other for enhanced visibility. The first 
spin fiip occurs if A* < 0, iV = 9 (a), and A* > 0, iV = 8 (c). 
No spin flip occurs if A* — (b). 



forms a half- filled shell with = 1, and the electron 
polarization P = 2/8, and A^ = 9 shows an open shell 
with Sz = 1/2, and P = 1/9. Because of d-shell over- 
turning, A^ = 10 forms a closed shell with Sz = 0, and 
P — 0/10. Here the energy difference between the KS 
eigenenergies of the electrons with spin up and down is 
zero, and A* « 0.6 meV (< Jem)- With increasing e- 
Mn couphng to Jem = 3.75 meV, corresponding to 5% 
doping, and within a range that spin polarization hap- 
pens only in d-shell, we find no magnetic transition for 
A^ = 8, since it was already fully polarized due to the 
spin Hund's rule, whereas for A^ = 9 and A^ = 10 we 
find transitions to s = 3/2 and s = 1 respectively. The 
evolution of quantized energies as a function of effective 
Zeeman energy and A* is shown in Fig. [4] For exam- 
ples A* < corresponds to d-levels of the 2D Gaussian 
potential with 7=1, whereas A* = (b), and A* > 
(c) correspond to d-levels in the 2D parabolic and 2D 
Gaussian potentials with 7 = (A* = A). 
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FIG. 5: The averaged magnetization per unit area (M^) as a 
function of number of electrons N a,tT = IK and Mn-density 



0.1 nm (a), and n„ 



0.025 nm"'^ (b) for non- 



n„ 

interacting (7 = 0, empty triangles) and interacting (7 = 1, 
filled triangles) electrons. At = 8, rim ~ 0.025 nm~'^, in 
2D Gaussian confining potential the ground state of the QD 
switches between ferromagnetic and antiferromagnetic states 
by e-e Coulomb interaction. 



We next turn to calculate the spatially-averaged Mn- 
magnetization per unit area A, (M^) = -j J p{Mz{p)) 
for the 2D parabolic and 2D Gaussian potentials as a 
function of N for zero (7 = 0), and full (7 = 1) e-e 
Coulomb interaction in Fig. [5l This calculation illus- 
trates how deviation from nonparabolicity of the confin- 
ing potential may lead to changes in the magnetic prop- 
erties of QDs. Many qualitative features indeed coincide 
for both of the confinements. 

A comparison in Fig.[5^a) for Um = 0.1 nm~^, confirms 
that there is similarity in magnetization for both Gaus- 
sian and parabolic confinements. Here Jem > A* (A) 
for 7 = 1 (0) and we consider all the electron numbers 
up to a full d-shell {N < 12). However, as we show 
in Fig. Efb) a nonparabolic confinement can introduce 
additional magnetic transitions in quantum dots. Here 
n-m = 0.025 nm^'^ corresponds to a magnetic doping of 



1.25 %. At N = 8 and 7 = in the 2D Gaussian poten- 
tial, we find that Jem < A and a transition to the AFM 
state. In contrast, at 7 = 1, we find Jem > A*, and thus 
the FM state is stable. 

Our findings illustrate the importance of Coulomb 
interactions and quantum confinement for the mag- 
netic ordering of carrier spin and magnetic impurities 
in (ll,Mn)VI quantum dots. We show that with elec- 
trostatic control of nonparabolicity of the confining po- 
tential, it is possible to control the magnetization of QDs 
even with a fixed number of electrons. While the choice of 
2D Gaussian confinement is a convenient form to param- 
eterize nonparabolic effects, other deviations from the 
parabolic confinement could also infiuence the magnetic 
ordering in quantum dots [28| . 
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